Optomechanical Entanglement in the Presence of Laser Phase Noise 



(N 

o 



m ■ 

Oh- 
^— > ! 

c : 

=5 

cr 



> 
oo 

00 

o 
o 



X 



R.Ghobadi 1 , 2 A.R.Bahrampour 2 and C.Simon 1 

' Institute for Quantum Information Science and Department of Physics and Astronomy, 
University of Calgary, Calgary T2N 1N4, Alberta, Canada 
2 Department of Physics, Sharif University of Technology, Tehran, Iran 
(Dated: January 27, 2013) 

We study the simplest optomechanical system in the presence of laser phase noise using the covariance matrix 
formalism. We show that the destructive effect of the phase noise is especially strong in the bistable regime. 
This explains why ground state cooling is still possible in the presence of phase noise, as it happens far away 
from the bistable regime. On the other hand, the optomechanical entanglement is strongly affected by phase 
noise. 
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I. INTRODUCTION 

Observing quantum effects for macroscopic objects like 
mechanical oscillators, if not impossible! 1]> 1S extremely hard 
due to the interaction of the system with its environment. Ac- 
cording to this explanation if one properly isolates the system 
of interest from its environment it would be possible to ob- 
serve bizarre quantum effects like nonlocality and entangle- 
ment at macroscopic level. One promising proposal to this 
end is to use an electromagnetic resonant system coupled to a 
mechanical oscillator. This type of system has been realized 
in various contexts and scales, including optomechanical sys- 
tems, electromechanical systems [2], Cooper-pair boxes [3], 
single electron transistors [4] and quantum dots [5]. Using 
this technique recently the first man made device has been 
prepared in the ground state and its coherent interaction with 
a qubit has been observed [6]. 

It has already been shown[7] that optomechanical systems 
can exhibit entanglement for small enough decoherence and 
strong enough optomechanical coupling. For typical optome- 
chanical systems cryogenic cooling is not sufficient to reach 
the quantum regime. To further cool the mechanical resonator 
one can use side band laser cooling by a red detuned laser. 
The laser serves both to cool the mechanical oscillator and to 
enhance the optomechanical coupling constant. In this way 
the strong coupling regime has been reached recently[8]. But 
the price to be paid is that a new source of noise is introduced 
in our system, namely laser phase noise (LPN). 

The conclusion of early studies of the effect of LPN jjj, 
based on modeling the LPN as white noise, suggested that 
with the current level of LPN ground state cooling is impossi- 
ble. This is in contrast with several experiments in which low 
phonon numbers have been observed. This discrepancy has 
been addressed in II 1 Ofl by considering a more realistic model 
for the LPN. 

Here we consider the effect of realistic LPN on optome- 
chanical entanglement. It turns out that optomechanical en- 
tanglement is very sensitive to this additional source of noise. 
It is particularly sensitive in the bistable regime, where the 
maximum entanglement is achieved in the absence of noise 

MM- 

The paper is organized as follows: Section II introduces the 
optomechanical system and establishes the notation. We ob- 



tain the governing equations in the presence of LPN in this 
section. In section III we introduce the LPN model. Sec- 
tion IV discusses the optomechanical cooling. We recover the 
well-known effect of LPN on the ground state cooling as a 
special case of our general results. Section V discusses the 
optomechanical entanglement in the presence of LPN. Sec- 
tion VI is a summary and conclusion. 



II. THE SYSTEM 

Our system is a high finesse Fabry-Perot cavity with one 
tiny end mirror. This mirror can move under the influence of 
the radiation pressure inside the cavity and in the same time 
undergoes brownian motion as a result of its interaction with 
the environment. The system is driven by a laser with fre- 
quency ll>l with power P. The Hamiltonian of the system is 
given by H — Hq + H mt where 



t t +- + tlLOm ,2 2\ 

Hq = n<D c a a H — — (q + p ) 



(1) 
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Eq.(l) describe the free Hamiltonian of the cavity mode and 
mechanical oscillator, respectively. u> c and a are frequency 
and annihilation operator of the cavity mode,respectively, a> m , 
q,p are frequency and dimensionless position and momen- 
tum operator of the mirror,respectively.The first term in Eq .(2) 
describes the optomechanical interaction in which Go = 

77 sjlmr ls tne sm gl e photon coupling constant. The second 
term in Eq.(2) describes the cavity pumping with Laser which 
its phase variation in time is given by (p(t) and E = Jj£*; 
where P, u>l are the input laser power and frequency respec- 
tively. 

In the frame rotating with frequency coit+ip{t) the equations 
of motion in the presence of damping and noise are 



q = u m p 



p = -u),„q - y m p + G a + a + £(t) 



(3) 



(4) 
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a — —{k + /Ao)fl + iaip + iG^aq + E + ^Jlica 



(5) 



where Ao = a> c - u>l and a m is the vacuum input noise. 
The mechanical and the optical input noise operators are 
fully characterized by their correlation function which in the 
Markovian approximation given by 



{a in (t)aUt')) = 6{t-t'). 



y m (2n + l)6(t - f). 



(6) 



(7) 



where n = [ejtp( jy ) — I] -1 is the mean thermal phonon num- 
ber and k B is Boltzmann's constant. The nonlinear Eqs. (3,4) 
can be linearized by expanding the operators around their 
steady state values (9, = (jJ + 60i(t) , where (9, - q ,p,a. 
The steady state solution of the system is given by q s = 



Gok,| 2 



p s = 0,a s = 



— 7, where a s ,q s ,p s are the values 
for cavity amplitude, position and momentum of mechanical 
oscillator,respectively. The last of these relations is a third 
order polynomial equation for a s , which can give rise to op- 
tomechanical bistabilitv lfl3i [l4ll . As can be seen from Fig. 1 
for strong enough input power the intracavity power admit two 
solutions. Fig. 1 shows the hysteresis loop for the intracavity 
power. Consider P cm , initially on the lower stable branch (I in 
Fig. 1). By increasing the input power we approach the end 
of the first stable branch. Increasing the input power beyond 
the end of the stable branch causes the power inside the cavity 
to switch to the second stable branch (II in Fig. 1). For in- 
put power larger than this value the cavity power is given by 
the upper branch. If the input power is decreased again, the 
cavity power follows a hysteresis loop. The parameter which 
describes the distance from the end of the stable branches is 
called bistability parameter which is defined bv lfT2l 



7=1- 



G A 



OJ m (_K 2 + A 2 ) 



(8) 



which is a real number between zero and one as long as the 
system has a stationary state[|7|]. As can be seen from Fig.l, rj 
decreases when approaching the bistable regime and becomes 
equal to zero at the end of each stable branch. 



Introducing amplitude and phase quadratures as X 



and Y = 



V2« 



and corresponding noises X,„and T„ 



6q = a) m 5p 



Sa+Sa* 
V2 

we have 



(9) 




FIG. 1: Bistability of the intracavity power with respect to the input 
power jT^]- The solid and dotted lines correspond to the stable and 
unstable branches respectively. The inset shows the bistability pa- 
rameter t] for the two stable branches. The end of each stable branch 
corresponds to rj = 0. 



where G = y/2Goa s , A = Ao - Goq s are the en- 
hanced optomechanical coupling rate and effective detun- 
ing. Introducing u T (t) = (q(t), p(t),X(t),Y(t)) and n T = 
(0,£(f), V2kX,„, ^2i<Yi n + V2a s £), Eq.(8-ll) can be written 
in a compact form 



u(t) = Au{t) + n(t) 



where 



A = 





-CD, 


G 



a) m ^ 
-7m G 



-K A 
-A -K J 



(13) 



(14) 



The solution of Eq.(13) is given by u(t) = M(t)u(Q) + 
J dsM(s)n(t - s) where M(t) = exp{Ai). The system reaches 
steady state if real part of the eigenvalues of A be negative 
which is the case if < tj < 10]. 

Since the initial state of the system is Gaussian and the dy- 
namical equation of system is linear in creation and annihila- 
tion operator both for cavity and mechanical mode the state of 
the system remain Gaussian at all time. A Gaussian state is 
fully characterized by its covariance matrix which is defined 
at any given time t by V, 7 (f) = ' 2 1 . Using the so- 
lution of Eq.(13) one obtain the following equation for steady 
state covariance matrix 



5p — -u> m 6q - j m 5p + G6X + £(t) 



(10) 



AV + VA' 



-D. 



(15) 



where the diffussion matrix is given by 



6X = -k6X + A6Y+ ^lKX in (11) 
6Y = -k8Y - A6X + G5q + V2^F,„ + ^2a s Cp (12) 



D tJ = ^f ds(M ik (s)({n k (s), bj(0)}>+ Jf^aXfaO)), n k (s)})). 

(16) 

Using Eq.(6,7), we get D = Diag[Q, y m (2n + 1), k, k + N] and 
N = L dsM44(s){(p(s)<p(0)) describes the contribution of the 
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LPN to the decoherence. Using Eq.(15) with the modified 
diffusion matrix the effect of arbitrary LPN on the behaviour 
of the optomechanical system can be obtained. 

III. THE NOISE MODEL 

The dynamics of phase for an ideal single-mode laser far 
above the threshold is given by itlrSjllTll 

ip + y c 4> = Z v {i) (17) 
in which ^(f) is a Gaussian random variable obeving |fl6l 



)) = 2y 1 c T L 5(t-t). 



(18) 



Here Ti and y^ 1 describe the laser linewidth and finite corre- 
lation time of phase noise. Using Eqs. (7) and (18), the phase 
noise spectrum is given by 



2T L 

1 + fi£ 
1 + rr 



(19) 



Note that for y c — > oo one recovers the white noise, while for 
finite correlation time and u > y c the power spectrum of LPN 
is greatly reduced. This means that a white noise model may 
greatly overestimate the effects of LPN. For this specific noise 
model we obtain 



N = 2a]y c r L 



f 

Jo 



dsMi,4(s)e 



(20) 



The general expression for N is very complicated and will 
not be included here. But it is still possible to get some in- 
sight about the behaviour of N qualitatively. Mi^is) is one 
of the elements of the matrix M(s) = exp(As). Far from the 
bistable regime, all eigenvalues of A are significantly smaller 
than zero, and all elements of M(s) decay with s, effectively 
limiting the range of the integral in Eq. (20). However, as one 
approaches the bistable regime, one of the eigenvalues of A 
approaches zero, and the range of the integral increases until 
it is limited only by the exponential e~ 7cS . One thus expects N 
to take its maximum values in the bistable regime. 

Fig. 2 shows N/k as a function of effective detuning. It can 
be seen from Eq.(20) and also in Fig. 2 that the noise is pro- 
portional to laser bandwidth, i.e F^. One also sees that N/k is 
maximum for 77 ~ 0, in agreement with the above argument. 
We would like to emphasize that the above argument is not re- 
stricted to the specific LPN model which we use here. In fact 
as long as the system admits a stationary solution (this guar- 
antees that M44(s) is a monotonic decreasing function of s), 
the above argument is valid for any LPN model with a finite 
correlation time. 



IV. OPTOMECHANICAL COOLING 

For completeness we study the optomechanical cooling in 
the presence of LPN. Solving Eq. (15) one can obtain the 
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FIG. 2: Plot of N/k and 77 as a function of effective detuning. We 
consider a Fabry-Perot cavity with length L = Imm and k = 1.4w m , 
driven by a laser with A = 81 Oram and input power 50m W. The 
mechanical oscillator frequency, damping rate and mass are 10 
MHz, 100 Hz and 5 ng respectively. The solid line corresponds to 77 
(multiplied by 60). The dashed and dotted dashed curves correspond 
to T L = 100 Hz and 30 Hz respectively. 



phonon number n m - 



V11+V22-1 



The quantum limit of phonon 



number can be derived by assuming a high mechanical quality 
factor and low temperature environment,!. e. y 1 » 1 and 
r^- >> 1. For rj ~ 1,k « a> m and A = cj m we find 

v 2 N 

(21) 



N 
+ — 
4k 



N = 



2a 2 J L y c (y c + k) 

w m + (y c + K ) 2 

Assuming a> m » y c » k we find 

2 2 
k a- 

" m = 4ZZ + TK S * i0Jm) 



(22) 



(23) 



which is identical to Eq. (23) in IToll if one add the effect of 
residual thermal occupation. From Eq.(19,23) its clear that as 
long as oj„, » y c the LPN effect on ground state cooling can 
be small. For y c « k the second term in Eq.(23) becomes 

a 1 

jj-SyicJm). Intuitively the possibility of ground state cooling 
in the presence of LPN can be understood as a result of the fact 
that ground state cooling happen around 77 ~ 1, far from the 
bistability regime, and thus in a region where the effect of LPN 
is small. See also the discussion of ground state cooling from 
the point of view of the bistability parameter in Ref . IU2I1 . 

V. OPTOMECHANICAL ENTANGLEMENT 

We use logarithmic negativity as the measure of entangle- 
ment, which is defined as [22] 



Em = max{0, -/n(2v m ,„)}. 



(24) 
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FIG. 3: Entanglement as a function of the bistability parameter 77 and 
the normalized detuning A/a>„, for three different values of the laser 
linewidth T L , from top to bottom: T L = 0, 10 Hz and 100 Hz. The 
other parameters are the same as for Fig. 2. 



where v,„,„ is the smallest symplectic eigenvalue of the 
partially transposed covariance matrix given by v„„„ = 
^ i-Vg-4<fo^ where 2 = detA + detB _ 2detC, and we repre- 
sent the covariance matrix in terms of 

V = [£t C b \ (25) 

Fig. 3 shows the entanglement as a function of the bistabil- 
ity parameter and the detuning for three different values of the 
laser linewidth. In the absence of phase noise, entanglement is 
maximal in the bistable regime, i.e. for 77 ~ 0. However, this 
changes dramatically as soon as the laser linewidth is non- 
zero. The entanglement goes to zero for 77 = for a linewidth 
as small as Fl — 10 Hz. For F^, = 100 Hz only a small amount 
of entanglement survives, and the region of maximum entan- 
glement is comparatively far from the line 77 = 0. It is no- 
ticeable that for the F^ = 100 Hz entanglement survives only 
for A/a> m > 2. This can be understood from Fig. 2, because 
the noise term N becomes very small for these values of the 
detuning. 



VI. SUMMARY AND CONCLUSION 

We studied a generic optomechanical system in the real- 
istic situation where the input laser bandwidth can not be 
ignored. We found that the LPN contribution to the deco- 
herence, characterized by N, is particularly significant in the 
bistable regime (77 ~ 0), and significantly suppressed else- 
where. This explains why ground state cooling is still possi- 
ble, as it happens for 77 ~ 1 . In contrast, optomechanical en- 
tanglement in the absence of LPN is maximal in the bistable 
regime. As a consequence, both the optimum region for the 
observation of entanglement and the amount of entanglement 
that can be achieved are strongly affected by LPN. 

Note added. When this work was completed, we became 
aware of a recent related paper ll23l . In comparison, the au- 
thors of that paper use a different noise model, and they treat 
the laser phase and amplitude as additional dynamical vari- 
ables. In contrast we keep the number of dynamical variables 
fixed and treat the LPN via an additional term in the diffusion 
matrix. 
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